Noise-enabled precision measurements of a DufRng nanomechanical resonator 
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We report quantitative experimental measurements of the nonlinear response of a radiofrequency 
mechanical resonator, with very high quality factor, driven by a large swept- frequency force. We 
directly measure the noise-free transition dynamics between the two basins of attraction that appear 
in the nonlinear regime, and find good agreement with those predicted by the one-dimensional 
Duffing equation of motion. We then measure the response of the transition rates to controlled levels 
of white noise, and extract the activation energy from each basin. The measurements of the noise- 
induced transitions allow us to obtain precise values for the critical frequencies, the natural resonance 
frequency, and the cubic nonlinear parameter in the Dufhng oscillator, with direct applications to 
high sensitivity parametric sensors based on these resonators. 






cd 



X3 

o 



> 

OO 
(N 

o 
o 



o 



X 



Doubly-clamped mechanical resonators have recently 
been the subject of much attention, due to the ability 
to make very high frequency, high quality factor res- 
onators, with applications in weak force and small mass 
detection, frequency stabilization, and possibly quantum 
computation^i^i^i'^i^i^i'''i^i^i^°i^^i^^i^^. The limit for para- 
metric sensing is often set by the precision with which 
a resonator parameter, such as the mass, can be mon- 
itored, limited typically by measurement and intrinsic 
noise sources. Here we show how one can use the in- 
trinsic nonlinear response of these resonators, and the 
addition of external broadband noise, to significantly im- 
prove the measurement precision of two such parameters, 
the resonance frequency and the cubic nonlinearity. This 
has direct implications for the ultimate sensitivity of such 
parametric sensors. 

At large drive amplitudes, doubly-clamped resonators 
exhibit a bistable response quantitatively similar to that 
of the Duffing oscillatorii^ii^. The motion in the fun- 
damental mode of a doubly-clamped beam is thus well- 
approximated by the Duffing equation, which for a natu- 
ral resonance frequency fig and quality factor Q, driven 
at frequency f2, has the form 
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where Y denotes the displacement amplitude of the mid- 
point of the beam, M denotes the mass of the beam, B 
the amplitude of the external driving force, and Bnoise(t) 
the stochastic forcing function due to thermal and ex- 
ternal noisrfiii^. This equation assumes that the beam 
oscillates in the mode with natural frequency Qq, that 
the displacement amplitude Y(t) is the only relevant de- 
gree of freedom, and that the equation of motion includes 
only the third-order nonlinearity, with strength K . 
The displacement Y{t) in Eq. (Q can be written as 



Y{t) = Ui{t) cos(fii) + C/2(i) sin(f7i), 



(2) 



in terms of the two quadrature amplitudes f/i,2(i)- For a 
high Q system driven at frequency 17 near fio, the slowly- 
varying envelope approximation can be usedi^ti^, where 



the functions f/i,2(i) are replaced by their slowly varying 
averages, ui,2(i)j respectively. 

In the absence of noise, the average functions ui^2{t) 
satisfy the equations of motion 
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The Duffing oscillator exhibits one stable state for 
small drive amplitudes B, while above a critical ampli- 
tude Be a bifurcation occurs, creating two stable basins 
of attraction. One basin corresponds to larger displace- 
ment amplitudes, and is stable for drive frequencies up 
to an upper critical frequency vu {v — f2/27r), deter- 
mined by the drive amplitude B. The other stable basin 
has smaller displacement amplitude, and is stable for fre- 
quencies down to a lower critical frequency ul, also deter- 
mined by the drive amplitude. The stable attractors are 
found by setting all time derivatives in Eq. (jS)) to zero 
and solving for mi,2, yielding three equilibrium points. 
Two of these equilibrium points are stable foci, and the 
third is a metastable separatrix. 

A transition between the two basins occurs in the ab- 
sence of noise when the energy barrier separating them is 
reduced to zero, by changing either the drive amplitude 
or the drive frequency. In the presence of noise, however, 
the Duffing oscillator will exhibit stochastic transitions 
between the two basins. For weak noise signals, the tran- 
sitions occur only near the critical frequencies vl.u^ while 
as the noise power is increased, the separation between 
the upper and lower transition frequencies is effectively 
reduced. 

Here we make detailed measurements of the nonlin- 
ear dynamics of a doubly-clamped beam, investigating 
both the dynamical motion and the change in the inter- 
basin transition rates due to broadband noise. Our ex- 
perimental system comprises a pair of doubly-clamped 
beams of single-crystal aluminum nitride, with dimen- 



sions 3 X 0.2 x 0.14 ^m'^, oriented perpendicular to one 
another and fabricated together on a chip of single-crystal 
Si. The fabrication technique is described elsewhereii. 
The chip was placed in the vacuum bore of an i? = 8 
T magnet at 4.2 K, with one beam (the active beam) 
oriented perpendicular to the field direction, the other 
(reference) beam parallel to the field. Magnetomotive 
actuation and displacement detection was used to drive 
the active beami, where the parallel orientation of the 
reference beam decouples it from the drive force (see Fig. 
^. The active beam had a natural resonance frequency 
1^0 = no/2TT = 92.9 MHz, a quality factor Q = 6750, and 
a critical drive power for inducing the hysteretic bifurca- 
tion of -61 dBm. Using the beam resistance of 11 fl, this 
corresponds to a critical drive force Be = 580 pN, and a 
midpoint displacement of 18 nm. 

Measurements were made with a radiofrequency 
bridge-", as shown in Fig. da). The rf drive signal is 
split by a 180° phase splitter, with the two phases passed 
through separate stainless coaxial cables of similar con- 
struction. The 180° phase-shifted signal is connected to 
one end of the reference beam, and the 0° signal con- 
nected to one end of the active beam. The other ends of 
the two beams are connected to a third coaxial cable that 
returns to room-temperature electronics. The bridge can 
be balanced in both amplitude and phase over the range 
of frequencies used in this experiment, and is typically 
tuned so that the electrical signal is proportional to the 
displacement-induced electromotive forcai. 

The signal measured in the experiment is the demod- 
ulated output of the bridge, giving the in-phase and out- 
of-phase quadrature signals I{t) and Q{t). These are 
proportional, to within a phase factor, to the average 
amplitudes ui^2{t)^^- In Fig. nib)-(d) we display the re- 
sponse of the active beam to a range of drive amplitudes, 
where the frequency is swept through the resonance for 
each drive amplitude; the hysteresis in the amplitude and 
phase response is in quantitative agreement with the re- 
sponse expected for a Duffing oscillator. 

In Fig. |2]we compare the measured quadrature ampli- 
tudes to numerical solutions of Eq. Q. In Fig. ^a) we 
show the calculated phase-space trajectories, and in (b) 
and (c) the experimentally measured trajectories. In Fig. 
I^fd) and (e) the time traces are shown for the switching 
transitions. The correspondence between the image in 
(a) and those measured in (b) and (c) is quite clear. It 
can be shown that as the drive frequency is varied, the 
stable points will follow a circle on the wi^2 plane, even 
in the presence of the Duffing nonlinearity. These circles 
are evident in Fig. ^d). 

We now turn to a discussion of the noise-induced 
transitions between the stable foci. The problem of 
thermally-activated escape from a potential landscape 
with a single basin of attraction is a thoroughly stud- 
ied probleroi^. The escape rate over a barrier of height 
Eb is given by F = a((3)z^o exp(— i^s/fesT), determined 
predominantly by the Arrhenius factor and less so by the 
Q-dependent prefactor a{Q). Our system differs from 
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FIG. 1: (a) Schematic showing circuit and active and refer- 
ence beams; dotted outline encloses cryogenic part of exper- 
iment. Box labelled 0, 180 is a 180° phase splitter, and that 
labelled A, allows adjustment of the amplitude and phase 
of the signal. The arrow indicates the orientation of the mag- 
netic field, (b) Hysteresis in amplitude versus drive frequency, 
for drive amplitudes from -68 to -53 dBm, in 1 dBm steps. 

(c) Hysteresis in phase for the same drive amplitudes as (b). 

(d) Hysteresis in u\ — U2 plane, plotted as / versus Q in di- 
mensionless units. 
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FIG. 2: (a) Numerically-generated phase-space flow for a 
drive force 9 dB above the critical point Be, and drive fre- 
quency 40 kHz above ^o/2tt. Flow begins near the separa- 
trix and evolves toward either focus, (b) Experimental phase 
space mean trajectory from focus 1 to focus 2 (8000 averages), 
(c) Data for phase space mean trajectory from focus 2 to fo- 
cus 1 (8000 averages), (d) and (e) Experimental time traces 
for the two switching transitions (8000 averages). 



this classic problem: Here, there is a basin of attraction 
about each of the two foci found on a Poincare map of 
the configuration space. Instead of a one-dimensional po- 
tential well, there is a quasipotential, with the dynamics 
governed by the noise energy at each point in the config- 
uration spacers. The equivalent activation energy, Ea, 
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FIG. 3: (a) Switching histograms h{v) for different noise pow- 
ers, with B —-56 dBm, 5 dB above the critical point, for tran- 
sitions from focus 1 to 2. (b) Transition rates r{i') extracted 
from the switching histograms, (c) Calculated energy barrier 
Ea{i^) extracted from the transition rates and the variation in 
noise power. The noise power was varied from -127 dBm/Hz 
to -113 dBm/Hz. 



for transitions between the foci, is found by integrating 
the minimuni available noise energy over the trajectory 
between the foci. 

Transitions were induced by using an external broad- 
band white noise signal, combined with the radiofre- 
quency drive signal using a directional coupler, to gen- 
erate a signal that included both the drive signal B and 
the noise signal S„. Typical noise powers ranged from 
-130 to -100 dBm/Hz. The drive signal itself was pro- 
duced by a source with very low phase noise; with no 
additional noise power, transitions were still induced by 
this remnant phase noise, to which the resonator is very 
sensitive. The thermal noise of the circuit, and the me- 
chanical noise associated with the finite resonator Q, are 
estimated to be 70 dB below the source phase noise, and 
were too small to induce measurable transitions in the 
system. 

Transition histograms were measured by applying a 
drive signal to the resonator above the critical value, 
preparing the resonator in one of the two basins of at- 
traction, and monitoring the switching transitions to the 
other basin. We measured histograms of the switching 
probability per unit time, h(t), by sweeping the drive 
frequency vit) = Q{t)/2Tr at a constant rate s = dv/dt, 
and recording the drive frequency at which a transi- 
tion occurred. This is a technique that has been ex- 
tensively used for measuring switching distributions in 
current-biased Josephson junctions^'^. The transition 
rate T{v) is extracted from the histogram h{t) using 

T{v{t))^{i-f_^h{t')dtr'sh{t). 

In Fig. ISJa) we display a set of histograms h{v{t)); 
higher noise powers shift the peak switching frequency 



and also broaden the distribution. In Fig. Olb) we show 
the transition rates extracted from these histograms, 
demonstrating the rapid increase in transition rate as 
the noise power is increased. We then extract the quasi- 
activation energy EAii'), by inverting the thermal activa- 
tion expression T{i>) = roexp(— _E^(z^)/fcBTe/^), where 
the effective temperature Te// is proportional to the noise 
power, and the prefactor Fq is related to the Kramers low- 
dissipation form-^^, Fq ~ vqIQ. We note that in this tech- 
nique, the histograms are only logarithmically sensitive 
to Fq, so that a precise determination is difficult. In Fig. 
ISfc) we display the activation energy Ea{v) extracted 
from the histograms, showing the expected decline in the 
barrier energy as the drive frequency approaches the crit- 
ical frequency. The distributions shown in Fig. Otb) are 
seen to collapse onto a single curve EAiy). In Fig. 0{a) 
we show a collection of experimentally measured EAiy) 
curves, measured for transitions from focus 1 to 2 and 
from 2 to 1, for different drive amplitudes. 

We calculated the activation energies numerically. The 
dynamic solutions to Eq. J^l without noise give the re- 
laxation from the separatrix to one of the foci. During 
a noise-induced transition, the system is excited from 
a basin near a focus towards the separatrix, which it 
crosses and then relaxes to the other focus. There is 
an infinite number of possible trajectories that allow a 
transition. Given a specific trajectory, it is possible to 
calculate the contribution of the noise force using Eq. 
(Q. The total energy transferred to the resonator for 
a particular trajectory is found by integrating the noise 
power along that trajectory, thus yielding the effective 
quasienergy between the foci. The energy transferred 
is thus an action-like quantity, and the most likely es- 
cape trajectory is that which requires the minimum ac- 
tion. The action-like integral S of the system is then 

The most likely path yo(^) minimizes the integral S. 
Because the separatrix is a saddle point, the extremal tra- 
jectory will most likely travel near the separatrix. The 
oscillator will naturally evolve from a point near the sepa- 
ratrix to either focus, without contributing to the action- 
like integral, as this relaxation does not require a noise 
term. Only when the oscillator is evolving against the 
dissipative flow field, from a focus toward the separatrix, 
will it contribute to the action integral. 

We used a numerical minimization of the possible tra- 
jectories Y{t), using S* as a test function to approach the 
extremum trajectory Yo{t). Minimum trajectories were 
calculated for different drive frequencies and amplitudes, 
yielding the energy barrier as a function of the drive am- 
plitude, shown in Fig. ^b). We find good agreement (to 
logarithmic accuracy) between the measured and calcu- 
lated energy barriers. 

Near the critical drive power Be, analytic forms in- 
dicate that the energy barrier should have a quadratic 
dependence on the offset from the critical frequency 
(v — Vc)^ (where Vc = i^u,l)^ ■ This quadratic depen- 
dence is shown in Fig. |3tc) for one drive power, and for a 
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FIG. 4: (a) Measured quasienergy barrier for a drive force 
ranging from 1 dB to 7 dB above the critical force Be- (b) 
Numerically-calculated energy barriers between foci 1 and 2, 
assuming a Duffing oscillator with the same frequency, quality 
factor, and nonlinearity as measured from Fig. (c) Log-log 
plot showing Ea oc {v — VcY dependence of the energy barrier 
near the critical point, (d) Vc versus drive amplitude B/B)c. 
At large amplitudes the data diverges from the analytic form. 



range of different drive powers in Fig. EJc). This allows 
a determination of Vc for a given drive power; our typical 
histograms yield an uncertainty of /!^Vc/vc ~ 3 x 10~^. 
By comparing the experimentally-observed dependence 
of Vc on drive power with that obtained from numerical 
analysis of the Duffing equation, shown in Fig. ^d), we 
can extract the natural resonance frequency vq and the 
coefficient of nonlinearity K] the former is the intercept 
of the curves shown in that figure, and the latter related 
to the slope of the curves at small drive powers. We find 
i/Q = 92887360 ±10 Hz and K = (3745 ±4) x 10" N/m^. 
The frequency measurement represents a relative preci- 
sion of /^vq/vq k, 1.1 X 10~^. This level of frequency 
resolution has significant implications for e.g. mass sens- 
ing with mechanical resonators^iSS. 

The measurements described above were made in the 
small-to- moderate noise limit, with noise energies much 
less than the energy barrier. At higher noise powers, 
the hysteresis due to the nonlinear response can actually 
be quenched, by rapid noise-induced transitions between 
the two foci. This quenching is demonstrated in Fig. 
As the noise power is increased, the area of the hystere- 
sis loop grows visibly smaller, until, at the highest noise 
powers, the switching is no longer hysteretic. In this 
limit, the oscillator generates random telegraph signals 
as it makes transitions from one focus to the other. The 
spectrum of the random telegraph signal is related to the 
transition rate of the oscillator. 
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FIG. 5: Amplitude hysteresis plots, for no noise power (bot- 
tom), with the drive amplitude set at -59 dBm, 2 dB above 
the critical point. The noise power was increased in 2 dB 
steps for each succeeding frame. At the largest noise power, 
the hysteresis is quenched. 



In conclusion, we have measured the configuration 
space trajectories, and the transition rates, between the 
bistable states of a nonlinear radiofrequency mechanical 
resonator. These measurements are in good agreement 
with numerical simulations based on the Duffing oscilla- 
tor equation of motion. Detailed analysis of the noise- 
induced switching transitions allows a quantitative mea- 
surements of the energy barrier between the stable foci, 
and provides a highly sensitive measurement of two key 
resonator parameters, the resonator natural frequency 
and the nonlinear parameter. 
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